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Certain inflation models undergo pre-heating, in which inflaton oscillations can drive parametric 
resonance instabilities. We discuss several phenomena stemming from such instabilities, especially 
in weak-scale models; generically, these involve energizing a resonant system so that it can evade 
tunneling by crossing barriers classically. One possibility is a spontaneous change of phase from a 
lower-energy vacuum state to one of higher energy, as exemplifled by an asymmetric double-well 
potential with different masses in each well. If the lower well is in resonance with oscillations 
of the potential, a system can be driven resonantly to the upper well and stay there (except for 
tunneling) if the upper well is not resonant. Another example occurs in hybrid inflation models 
where the Higgs fleld is resonant; the Ifiggs oscillations can be transferred to electroweak (EW) 
gauge potentials, leading to rapid transitions over sphaleron barriers and consequent B-l-L violation. 
Given an appropriate CP-violating seed, we find that preheating can drive a time- varying condensate 
of Chern-Simons number over large spatial scales; this condensate evolves by oscillation as well as 
decay into modes with shorter spatial gradients, eventually ending up as a condensate of sphalerons. 
We study these examples numerically and to some extent analytically. The emphasis in the present 
paper is on the generic mechanisms, and not on specific preheating models; these will be discussed 
in a later paper. 

PACS numbers: 98.80.Cq, 98.8Q.-k, 11.15.-q , 05.70.Fh 



I. INTRODUCTION 

There are well-known reasons to believe that inflation 
took place and was followed by reheating to some temper- 
ature r^. Before a thermal equilibrium was reached, the 
coherent oscillations of the inflaton could create the envi- 
ronment in which a resonant non-thermal production of 
particles could rapidly transfer energy from the inflaton 
to the other fields. This stage, known as preheating [Q, 
has been a subject of intense studies. In particular, it was 
argued that both non-thermal phase transitions |^] and 
the generation of baryon asymmetry |^|^ could occur 
during preheating. 

We will describe two new field-theoretical phenomena 
that can be caused by coherent oscillations of the infla- 
ton. One is a new example of a phase transition driven by 
the coherent oscillations of the inflaton. This transition 
has an unusual feature that it can start in a lower-energy 
ground state and end in a higher-energy metastable vac- 
uum. We discuss this in Section II. 

In Section III we describe resonant generation of a 
fermion density through anomalous gauge interactions 
that can be the basis for baryogenesis. In contrast 
with the earlier work, where the analyses were based 
on analogies with thermal sphalerons j^J^ or topologi- 
cal defects Q , we construct an explicit solution that can 
be though of as a condensate of sphalerons. We show 



that the evolution of this solution can lead to a resonant 
growth of Chern-Simons number density. 



II. PHASE TRANSITIONS AT PREHEATING 

The properties of the physical vacuum and the par- 
ticle content of the universe are determined by physical 
processes that took place in a hot primordial plasma. 
Theories of particle interactions beyond the Standard 
Model allow for different types of physical vacua. For 
example, an SU(5) Grand Unified Theory (GUT) al- 
lows three possibilities for the ground state, in which the 
gauge group that remain unbroken is, respectively, SU(5), 
SU(4)xU(l), or SU(3)xSU(2)xU(l). If low-energy su- 
persymmetry is assumed (to assure the gauge coupling 
unification and to stabilize the hierarchy of scales), these 
three ground states are degenerate in energy up to small 
supersymmetry breaking terms ~ TeV. Therefore, any of 
these potential minima could equally well be the present 
physical vacuum. The evolution of the universe shortly 
after the Big Bang must have chosen SU(3)xSU(2)xU(l) 
vacuum over the others. The phenomenon we will discuss 
can provide a new solution to the old puzzle related to 
breaking of a SUSY GUT gauge group. The same pro- 
cess can have important consequences in other models 
with several competing (metastable) vacua, for example. 
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in the minimal supersymmetric extention of the Standard 
Model (MSSM). 

Let us consider an inflaton (f> interacting with a "Higgs 
field" X through a coupling of the form A$^x^X or 
lJL<^X^Xi or both. Let us assume that the effective po- 
tential V{x, *&) has two non-degenerate minima, for ex- 
ample at (x) — ±w, ($) = Wj, and that the mass of 
the X particle is not the same in both minima, that is 
d^V{v,v,)/dx^ ^d'^V{-v,v,)/dx^- 

At the end of inflation, the system can occupy the 
lowest-energy state with (x) — ~v. During preheat- 
ing, the inflaton oscillates around its VEV, = 
Vj + $0 cosuit. In general, this induces a time-dependent 
mass for the Higgs field x through the couplings and 
A. The equation of motion for the homogenous (zero- 
momentum) mode of the field x is 

d 

X + 3iIx+ ^^(x^-f. +*ocoswt) = 0, (1) 

where H is the Hubble constant.^] In Fig. |^ and Fig. || 
we show two examples of time-dependent effective poten- 
tials. 

The potential V{x,Ht)) = (x^ - v'^fn + 

0. 4cos5.6wi] + 0.1ux(3w^ — x^) depicted in Fig. |l] has 
two classical solutions, x = ~^ and x = ^- Naively one 
could expect that the lowest-energy solution x = —v cor- 
responds to the vacuum state. This is not necessarily 
the case, however. Since the mass of the x field is time- 
dependent, the solution x(i) = may be unstable with 
respect to small perturbations. At the same time, the 
other solution, x(i) = +v may be stable. If this is the 
case, the classical system is attracted to the trajectory 
X(<) = +v. 

In the vicinity of the global minimum, for | (x+^') lv\ ^ 

1 , the equation of motion (|l|) is a Mathieu equation that 
has rapidly growing solutions for some values of lj, 
and m*^~) = d'^V{—v,Vj)/dx^- The inflaton frequency 
changes with time and can enter in resonance, at which 
point (Ix(^)l begins to grow exponentially. This kind 
of solution of equation (1), with H=Q and the potential 
of Fig. |l|, is shown in Fig. |[ At some point it crosses the 
barrier and begins oscillations around a different poten- 
tial minimum, (x) = ^v. However, the mass of the x par- 
ticle near (x) = -\-v is m^^\ different from m!'~\ There- 
fore, the system may go out of resonance after crossing 
the barrier. There are no growing solutions in the vicin- 
ity of the seconds minimum, and the oscillations die out 
with (x) = 

If the tunneling rate between (x) = -\-v and (x) = —v 
is negligible, the classical evolution shown in Fig. || de- 
scribes a phase transition into a metastable false vacuum. 



^In weak-scale preheating the Hubble constant is negligiblly 
small. For GUT-scale preheating it is not, and it could play 
an important role in helping to scan resonant bands. 
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FIG. 1. The time-dependent potential 

y"(x,$(t)) = (x^-w^)^[l + 0.4cos5.6i;t]-h0.1tJx(3tJ^-X^) that 
has two non-degenerate minima and a time-dependent barrier 
height. The masses of the x particles are also time-dependent 
and are different in the two minima. 
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X/i; 

FIG. 2. Another time-dependent potential. The heights 
of the two vacua oscillate. 

This example shows that the ground state at the end 
of inflation does not necesarily correspond to the global 
minimum of the potential. Instead, during the preheat- 
ing, a false vacuum can be populated if the true vacuum 
entered in resonace while the false vacuum did not. 

Both Grand Unified Theories and supersymmetric ex- 
tentions of the Standard Model predict the existence of 
local minima in the effective potential. The tunneling 
rate between these minima can be extremely low and 
their lifetimes can easily exceed the present age of the 
universe. For example, the effective potential of the 
MSSM can have a broken color SU(3) in its global min- 
imum, while the standard, color and charge conserving 
vacuum is metastable. For natural and experimentally 
allowed values of the MSSM parameters, the lifetime of 
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FIG. 3. Classical solution of the equation of motion in the 
potential of Fig. |^. The evolution begins near the unstable 
classical trajectory = ~v and is driven towards a stable 
classical solution x{i) = +'"■ In quantum theory, if the tun- 
neling rate between the two vacua is small, a phase transition 
to a metastable vacuum takes place. 



this false vacuum can be much greater than lO^*' years . 
If the reheat temperature after inflation was not much 
higher than the electroweak scale, this metastable mini- 
mum could be populated in the way we have described. 

Breaking a SUSY GUT gauge group and choosing be- 
tween the nearly degenerate minima is problematic in 
non-inflationary cosmology Let us consider a SUSY 
SU(5) GUT for example. The minima with unbroken 
SU(5), SU(4)xU(l), and SU(3)xSU(2)xU(l) groups are 
nearly degenerate, split only by supersymmetry breaking 
terms of the order of a TeV. Why did the universe end 
up in the vacuum with the lowest symmetry? 

Finite temperature corrections (if relevant, which may 
not be the case for preheating) make the SU(5) minimum 
lowest in energy because it has a higher number of degres 
of freedom. The subsequent thermal evolution of the 
potential makes tunneling into a Standard Model vacuum 
impossible |8| even if it becomes the global minimum at 
temperatures below 1 TeV. Supergravity splits the three 
minima by a negligible amount and in such a way that 
cosmological constant can by fine-tuned to zero only in 
the minimum with the higher energy while the other two 
minima have negative energy density Some of the 
proposed solutions |^ rely on assumptions about a strong 
gauge dynamics that seem somewhat implausible. 

If, however, inflation took place, the SUSY GUT vac- 
uum could be chosen in a phase transition of the kind we 
described. This appears to resolve a long-standing prob- 
lem concerning the breaking of the SUSY GUT gauge 
group. 



III. B+L VIOLATION 

As discussed in the Introduction, preheating oscilla- 
tions of the Higgs VEV can lead to two effects of interest 
for B-l-L violation. The flrst is that the sphaleron bar- 
rier itself oscillates, leading in principle to exponentially- 
sensitive oscillations of the sphaleron rate. The second, 
which we take up here, is that Higgs oscillations can res- 
onantly drive classical transitions over the barrier. 

Given an appropriate CP-violating seed, there are 
three stages to this classical resonant driving. In the 
first stage, the seed (which can be a source term or ini- 
tial conditions on the EW gauge potentials) drives large- 
scale generation of Chern-Simons (CS) number (topolog- 
ical charge) over spatial scales so large that spatial vari- 
ation can be ignored and only temporal variation saved 
in the classical equations of motion. In the second stage, 
gradients on shorter scales emerge, as a result of unstable 
growth of spatially-dependent perturbations. The seeds 
for these spatial modes might emerge from spinodal de- 
composition during inflation • As expected on gen- 
eral grounds from earlier preheating studies, the fastest- 
growing modes are those with large spatial scales. The 
third stage involves the generation of sphalerons, with 
spatial scales at the standard W-boson mass Mw- 

In all stages, we will ignore various back-reaction ef- 
fects; the expansion of the universe (in any case, neglible 
for weak-scale inflation); and damping produced by per- 
turbative decays (one order of aw higher than terms we 
keep). 

We discuss the first stage, which has important 
non-linear effects stemming from gauge-potential self- 
coupulings, both analytically and numerically. A partic- 
ular ansatz is used for the gauge potential, having only a 
time dependence. (This ansatz has been used some time 
ago [|l2| in a rather different scenario.) The analysis is in 
the same spirit as the conventional approach to low-order 
resonances in the Mathieu equation (see, e.g., Ref. @). 
But the lowest-order resonant-mode equations, two first- 
order differential equations, have a cubic non-linearity. 
Surprisingly, these coupled non-linear equations can be 
solved exactly in terms of elliptic functions. The non- 
linear terms not only provide a quartic potential opposing 
the growth of CS number but, as the CS number grows, 
the non-linear term also grows and drives the system off 
resonance. In effect, the cubic non-linearity causes the 
W-boson mass to increase. Interestingly, this increase 
can be offset by a secular increase of the frequency of 
Higgs oscillations, allowing resonance to be maintained 
for long periods of time with consequent large growth of 
CS number. 

In the second stage we include linear perturbations to 
the spatially-homogeneous equations of the first stage; 
these perturbations are considered to lowest order in spa- 
tial gradients, as characterized by a spatial momentum k. 
It is not possible to do a conventional dispersion-relation 
analysis of these equations, which have time-dependent 
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coefficients as determined by the temporal growth of the 
first-stage gauge potentials. We perform a numerical 
analysis of the three coupled linear differential equations 
which result. 

The third stage, in which gradients evolve to spatial 
scales ^ appropriate for sphalerons, is the hard- 

est to analyze, since an adequate treatment involves the 
solution of coupled partial differential equations with 
time-dependent coefficients. So we restrict ourselves to 
a crude, simple ffi'st step, reducing these partial differen- 
tial equations by a non-linear ordinary differential equa- 
tion for an approximate sphaleron-like mode. The rel- 
evant gauge-potential ansatz, first introduced by Bitar 
and Chang , was later used to analyze sphalerons 
above the EW phase transition, and was shown to have 
an effective barrier potential for the sphaleron which was 
numerically very close to that of a simple pendulum. We 
introduce an oscillating Higgs field, which causes this 
pendulum to be parametrically-driven. The ansatz is 
too simple to be used for anything more than estimating 
the rate of change of topological charge as the pendu- 
lum goes over its barrier once; we do this numerically. In 
principle, more complicated forms, representing multiple 
sphalerons, could be used, such as the ADHM construc- 
tion or those of 't Hooft or of Jackiw, Nohl, and Rebbi | |T^ ] 
multi-instanton form, suitably modified for Minkowski- 
space dynamics, but these have not yielded any insights 
for us. 

At all stages, the energy density associated with gen- 
eration of CS number is of order Airm'^ / g'^ , as would be 
appropriate for a gas of sphalerons with density ~ m^. 

A. First stage: homogeneous CS parametric 
resonance 

In what follows we always consider the Higgs field to 
have a given VEV, as determined by preheating effects. 
Introduce the conventional anti-hermitean gauge poten- 
tial, with coupling g included, by: 

gA, = (2) 

Our spatially-homogeneous ansatz is: 

gAo = 0; gA, - ( J)</,(0 (3) 

in which the group index is tied to the spatial index. By 
the conventional rules of charge conjugation and parity 
for the gauge potential, is C even, P odd, CP odd. 

It is important to note that this ansatz does not corre- 
spond to a non-vanishing VEV for an EW field. Gauge 
invariance alone is enough to ensure that there can be 
no expectation value coupling the space-time indices to 
group indices. 

One readily calculates the EW electric and magnetic 
fields: 



gE, = Go. - ( 9B^ = ^e»,/cG,fe = {^)cf>\ (4) 

Then one calculates the density W of Chern-Simons num- 
ber as: 



It is straightforward to check that W is the topological 
charge density Q, related to B-l-L violation through the 
anomaly equation. 

With the assumption of a given Higgs VEV, the equa- 
tions of motion for the gauge potential are: 

[D'', G^.] + M^{t){A, + {d,U)U-^) - 0. (6) 

Here the unitary matrix U represents the Goldstone 
(phase) part of the Higgs field. The mass term will be 
assumed to have the form: 

= 7T?{l + ecos{ujt)) (7) 

where m is the value of My/ with no oscillations. Later 
we will have occasion to consider a time-dependent fre- 
quency oj, but for now think of it as a constant. 

There must be some sort of CP-violating seed to pro- 
duce non-zero solutions of the equations of motion; these 
might stem from (long-scale) spatial gradients in the ma- 
trix C/, which acts as a source in equation (6), or from 
initial values of (j). Because the equations are unstable, 
there is little practical difference, and we choose to drop 
the U terms in the equations of motion, and then provid- 
ing a seed through initial values. Then there is a single 
equation for <j): 

ij) + 2cj)^ + {l + ecosrt)cj)^G (8) 

We have non-dimensionalized the equations of motion by 
measuring (j) in units of m and time t in units of . 
The parameter p has the value ujjra. 

Without the cubic non-linearity, this would be a stan- 
dard Mathieu equation. In the Appendix we analyze the 
coupled non-linear mode equations which arise for the 
lowest resonance (r = 2), and find that they can be solved 
exactly in terms of elliptic integrals. The qualitative fea- 
tures of this analysis are easy to anticipate: Equation 
(8) describes the motion of a particle in a quartic poten- 
tial. The oscillating term drives the particle up the wall, 
but eventually the particle gets out of resonance and falls 
back. This process can repeat quasi-periodically. 

We now turn to numerical analysis. Only a couple of 
examples will be reported, without attempting to choose 
parameters to correspond to realistic preheating scenar- 
ios. Parameters are chosen to illustrate specific effects; 
other parameter sets may show no interesting behavior 
at all. The runs reported here have initial values 

0(0) = 0.001; ^(0) = 0, (9) 
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FIG. 4. Time dependence of cf) for initial values as in 
equation (9) plus r = 2.3, e = 0.9. 

and large values of e, in the range 0.5-0.9. Because 
the equations are both non-linear and unstable, the fi- 
nal results are largely independent of the initial condi- 
tions as long as they are non-zero. As the initial val- 
ues are reduced, the time of onset of instability is some- 
times lengthened. Generally, there are two regimes (for 
constant frequency lo): The resonant regime, in which 
(j) grows to 0(1), and the non-resonant regime where 
stays small. We will only show the near-resonant cases 
in the figures. There is another regime in which uj grows 
secularly with time, and which leads to larger values of 
</>. 

Fig. ^ is a typical example of the behavior when u) 
or r is constant and fairly near resonance (in this case, 
r=2.2>). One sees that the envelop of grows to order 
unity, but periodically passes through zero and repeats. 
This is because e is near unity, and so system frequencies 
vary quite a bit, from 1 -I- e to 1 — e. 

Fig. 1^ shows the behavior when the frequency grows 
secularly. The onset of rapid growth is delayed because 
the system is originally fairly far from resonance, but 
then the envelop of |0| grows essentially linearly, coupled 
to the frequency change. The system is able to stay in 
resonance as (j) grows linearly, because the effective mass 
M of the (j) field (see the Appendix) is ~ -1-3(0^), 
and the effective ratio r = lu/M stays roughly constant 
if M grows at the same rate as u>. 

Fig. ^ shows the CS density (/j'^/Stt^ corresponding to 
the parameters of Fig. |[ The CS density grows roughly 
as , with (j) growing linearly in time as does to. 

With dimensionalized values of |(/)| ~ to, the CS num- 
ber density is of order 0.01 to^, corresponding to a large 
B-l-L density. Whether any of this CS density survives 
preheating to the reheating phase depends on whether 
there is a "graceful exit" to preheating generation of 
CS number, and this depends on factors not considered 
in this paper, such as back reaction, growth of finite- 
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FIG. 5. Behavior of (j) with the initial conditions of equa- 
tion (9), with b = 0.52, and a secularly-growing frequency 
r{t) = 1.7 + O.OOlt. 
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FIG. 6. The Chern-Simons density (see equation (5)) for 
the parameters of Fig. 

momentum modes, and linear damping by decay of the 
W-boson condensate. Additionally, there may be many 
domains large compared to but small compared to 
the Hubble size in which the values of 4> are uncorre- 
lated. This will reduce the effective global CS density by 
a factor of iV^/^, where N is the number of such domains. 
The ultimate fate of the processes considered here will be 
taken up in a future work, in which specific weak-scale 
preheating scenarios will be taken up. 



B. Second stage: evolution of spatially-varying 
modes 

Ultimately, there will be some CP- violating seeds with 
finite spatial gradients. Assuming that these seeds are 
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smaller than those for (p (as is reasonable following in- 
flation), these seeds will be driven by the time varia- 
tion of (/) as well as of the Higgs VEV. We will be con- 
cerned here only with the linearized equations for the 
spatially- varying modes, which we characterize in mo- 
mentum space. As is usual in preheating phenomena, 
the modes with the longest spatial scales (smallest k) 
grow fastest. 



The total vector potential is written as 



with 



Af^ taken from equation (3). The most general vector 

potential ga^ depending on a single vector k has time 
component 



IT ■ k 
2i 



(10) 



and space components 



1 



gaj = ^[('^j ~ ^3^ • + i^]abTakbP2 + kjf ■ fc/Sa]. 

(11) 

In equations (10,11) the hat indicates a unit vector, and 
ao,f3i are real functions of k^ and t. As before, we non- 
dimensionalize by dividing these functions by m, replac- 
ing t by mt, and k by k/m. Presumably the Fourier 
transforms in (10,11) vanish at an appropriate rate as 
/c ^ so as to change k into fc, although this will not 
matter in what follows. 

It is straightforward if lengthy to write out the lin- 
earized version of equation (6) (without the U terms): 



1 



Q 

= k^ 



2</>2 



1 + ecosri: 



Pi + Qf3i - 2k4>P2 + 2(/3i + /33 



= 0; 



(12) 



(13) 



P2 + Qp2 - 2/c0/3i + (j){ao - kPs) + 2(/)ao = 0; (14) 



f33 + Qf33 + k{ao - fc/33 - 20/?2) + 4/?i(/.2 = 0. (15) 

Even though these are linear equations for the modes 
ao,f3j they are impossible to solve analytically, because 
(j> is not an analytically-known function. We have solved 
them numerically, with various interesting results. Per- 
haps the most interesting is that these mode functions 
remain small and well-behaved for a long time, and then 
when (j) is large enough (of order unity) they show vi- 
olently unstable behavior. This is especially so for the 
case when the frequency uj is growing with time, as for 
Fig. ^. This is illustrated in Fig. 0, showing the evolu- 
tion with time of the linear modes for the parameters of 
Fig. ^ The mode functions were begun with initial val- 
ues which are 0.1 times those of (j) (see equation (9)). Of 
course, any other initial values can be gotten by scaling, 



since the equations are linear. The point is that when, 
for a given set of initial values of ao,f3j, these functions 
rise to be of 0(1), the whole problem becomes non-linear 
and presumably enters something like the third stage dis- 
cussed below. Note in Fig. ^ that the threshhold for non- 
linearity, with the given initial conditions, occurs at a 
(dimensionless) time of 200, which gives enough time 
to get big enough to be interesting (see Fig. ||). 



C. Third stage: sphalerons 

Eventually, momentum modes with A: ~ 1 will become 
prominent, and the condensate of CS number becomes a 
condensate of sphalerons. It is much more difhcult to de- 
scribe this stage, and we will only take a simple first step. 
This step consists of a drastic simplification of the kine- 
matics of a sphaleron coupled to a time-dependent Higgs 
field, reducing the dynamics to a single function X{t) 
as in Refs. [|l^,^. Write the most general spherically- 
symmetric gauge potential and Higgs phase matrix U in 
the form: 

1^ 1 
U ^ e^p[—f ■ f], gA„ = —f ■ TH2; (16) 
2 2i 



gAi = —[tiakTafki 

2ir 



1) - (Ti - fif ■ r)02 + nf ■ tHi. 

(17) 



The functions 

values of the angle 7 are zero at r 
We parametrize these functions as 



depend only on r, t. The asymptotic 
and TT at r = (X). 



Hi = 



2A 



A2 



H2 = - 



2rA 



(18) 



= 1 - 



A2 + r2 + a2 ' 
2r2 



A2 + r2 + a2 ' 
2rA 

'^^ " " A2 + r2 -f a2 • 

For details on the parametrization of 7 see For 
the present purpose one can just think of 7 as always 
equal to tt. In this parametrization the constant a is 
a size parameter (like that of an instanton) and A, the 
sole dynamic degree of freedom, depends on t. Generally, 
A is an odd function of vanishing along with its first 
derivative at t = 0. 

The electric and magnetic fields are: 



SB, = (S) 



4a2A 


A2 


+ T'^ + 




4a2 


A2 


_|_ J.2 _|_ 



(19) 
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FIG. 7. The evolution of the hnear modes for the param- 
eters of Fig. |5| Note the truncation on the vertical axis; at 
t ~ 200 

all the amplitudes are larger than one in magnitude. 



Note that these have the same space and internal sym- 
metry index dependence as does the 4> ansatz of equation 
(3). It is therefore natural to suppose that the (j) fields 



will transform (through the growth of spatial modes) into 
a condensate of sphalerons. Of course, in this condensate 
each sphaleron will be a translate in space and in time 
of the sphaleron exhibited here, which is centered at the 
space-time origin. 

With boundary conditions 



A(< = — oo) = —oo; X{t = +oo) = -|-oo 



(20) 



one readily verifies that, no matter what the dynamics 
of A as long as it is single- valued, the (Minkowskian) 
topological charge 



Q = / d'^xTrE ■ B 



(21) 



has the value 1. Indeed, if we replace A by t we get 
exactly the usual Euclidean one-instanton expression, 
which however is now being interpreted as a Minkowskian 
construct. 

The size coordinate a is not arbitrary, as it is for in- 
stantons in gauge theories with no Higgs field. As shown 
in if one goes to t = and sets A, A = there, 
the resulting ansatz in equations (16,17) is an excellent 
trial wave function for the sphaleron. Minimizing the 
Hamiltonian (for time-independent Higgs VEV) yields 
a — ^/i/2Mw) and a sphaleron mass Ms only a frac- 
tion of a percent higher than the true value, determined 
numerically, of 



Ms = 5.41 ^ . 

9 



(22) 



When the mass M-^^ depends on time, as in equation 
(7), we will continue to use the above value for a. It then 
happens that the parameters of the Hamiltonian depend 
on time (see [|l^ for the Hamiltonian as a function of 
a. A, A). 

As is further shown in p5| , one can trade the function A 
for a topological charge Q{t) defined by demanding that 
the kinetic energy term in the Hamiltonian is of the form 
(l/2)/Q^ with / independent of Q. The normalization 



A 



g = 0; A = +00 ; Q = 27r 



(23) 



makes the topological charge an angular variable. Nu- 
merical work shows that the potential energy is very 
nearly that of a pendulum, and that / — S^Ms/m? for 
some numerical constant ^. The resulting approximate 
Hamiltonian has the form: 



H ^ Ms[^Q^ - cosQ] 



(24) 



which has, as it must, the value Ms when (5 = 0, Q = n. 

Next one replaces M\y by its time-dependent value, 
as in equation (7). We have numerically investigated 
such driven pendulum equations. They lead to multiple 
transitions over the barrier, but we will not display such 
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solutions here. One reason is that the ansatz we use here 
is strictly tied to a unit change of topological charge, 
so that all that counts is the rate of making a single 
transition over the barrier. Just as for all the classical 
barrier-hopping solutions presented for the (j) ansatz^ the 
rate is 0{ijj), very much different from the timneling rate. 
(The tunneling^ate is also changed as the sphaleron mass 
oscillates; see [^.) 

To go further than this for a condensate of 
real sphalerons is extraordinarily complicated; each 
sphaleron, like the instanton to which it corresponds, 
has numerous degrees of freedom. Even if we restrict 
this to one degree of freedom (corresponding to A) for 
each sphaleron, it is not clear how to proceed. Nor is it 
clear how to modify known multi-instanton ansdtze such 
as ADHM or that of 't Hooft or Jackiw, Nohl, and Rebbi 
to express the real-time sphaleron dynamics in the 
presence of an oscillating Higgs field. 

IV. CONCLUSIONS 

In this work we have investigated two new mecha- 
nisms driven by preheating oscillations of, e.g., the Higgs 
field in hybrid inflation. The flrst mechanism, resonant 
barrier-crossing from a lower minimum to a higher min- 
imum (where there is no longer resonance), may explain 
some puzzles associated with the symmetry-breaking pat- 
terns of GUTs. This kind of transition could also popu- 
late a metastable SU(3) x SU(2) xU(l) vacuum in a super- 
symmetric extension of the Standard Model even if the 
global minimum of the potential breaks charge and color. 
(In the case of the MSSM, this posibility has direct inpli- 
cations for collider experiments 0.) The second mecha- 
nism, resonant barrier-crossing associated with B-l-L vi- 
olation, may lead to a condensate of sphalerons on time 
scales short compared to tunneling rates. Both effects 
require resonance with preheating oscillations to be ef- 
fective. We have not tried to construct "realistic" appli- 
cations of these mechanisms to specific preheating scenar- 
ios. We note, however, that in many cosmological mod- 
els, even if the initial conditions are far from resonance, 
the system evolves and reaches the resonance eventually, 
thanks to a change in the relevant parameters [Q. Such 
evolution is facilitated by either non-quadratic inflaton 
potential that causes a variation in the inflaton frquency, 
of by expansion of the universe and the associated Hub- 
ble damping (for GUT, not weak-scale preheating), or 
some other effects that can slowly drive a system into a 
resonance band. We leave the building of realistic cos- 
mological models for future work. 

Aside from such applications, there is still a good deal 
of work to be done to clarify these mechanisms. In the 
case of B+L violation, one can raise the following issues: 

1. How do the three stages (spatially- homogeneous 
potential, linear momentum-mode perturbations, 
sphaleron condensate) of Section HI evolve from 



the flrst to the last? This can only be answered by 
numerical work more extensive than we have yet 
done. 

2. The large-scale EW CS density we propose will 
have a projection onto Maxwell magnetic flelds 
carrying helicity (another term for Chern-Simons 
number). The spatially-homogeneous nature of 
these fields makes them quite different from ear- 
lier proposals (see and references therein) 
involving generation of Maxwell flelds in a ther- 
mal environment, with unacceptably small scale 
lengths to correspond to the scale lengths of 
present-day galactic magnetic flelds. Given suffi- 
cient inverse cascading of the nearly-homogeneous 
Maxwell flelds following from our preheating mech- 
anism (at EW time these fields must be limited in 
extent by the Hubble size), shows that EW- 
time Maxwell fields could indeed be the seeds for 
presently-observed galactic fields. We intend to in- 
vestigate this further. 

3. Can one make use of multi-instanton ansdtze such 
as those of ADHM, 't Hooft, or Jackiw, Nohl, and 
Rebbi to extend the Bitar-Chang construc- 
tion we have exploited in Section IIIC in order to 
understand quasi-analytically the formation of a 
sphaleron condensate? 

4. Are there (necessarily spin-dependent) quasi- 
resonant phenomena for the production of W- 
bosons by an oscillating Higgs field which are in 
any sense analogous to the very sharp resonant 
phenomena found by Cornwall and Tiktopoulos 
[ p9[ for spin- 1/2 charged particles in specific time- 
dependent electric flelds? 

To clarify this last point, Ref. [|9| found that it is pos- 
sible to have highly-resonant e^e^ pair production in a 
classical time- varying electric fleld of the proper time de- 
pendence. The sharply-resonant nature of the process 
can only happen for fermions, but in any case spin ef- 
fects, which might be available with gauge bosons, are 
important in overcoming the typical exp(— l/a) rate of 
pair production in classical flelds. 
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APPENDIX: ANALYSIS OF MODE EQUATIONS 

FOR (t> 

Here we give the analysis of the Mathieu-like but non- 
linear modal equations of Section HI. Just as for the 
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Mathieu equation, we write the non-dimensionalized 
in the form 



cj) = a{t) cos(rt/2) + b{t) sin(rt/2), 



(Al) 



(where, as in the main text, r = oj/m), leaving out all 
terms with higher frequencies. One verifies that the time 
dependence of a, 6 is 0(e), so that we can ignore second 
derivatives of these quantities. However, we will save the 
cubic non-linearities. 

Using equation (Al) in the equation of motion (8), 
saving only terms varying as cos(rt/2) and sin(rt/2), and 
dropping second derivatives yields: 



rb + a e 

U 2 



1 



(A2) 



(A3) 



To make contact with the linear Mathieu equation, let us 
temporarily replace the terms (3/2)(a^+fe^) by constants, 
and define an effective (non-dimensional, that is, scaled 
by m) mass M by: 



M2 = l + ^(a2+52)_ 



(A4) 



Assuming exponential growth, with a,h ^ exp(fit), gives: 



2r' 



(r2-4M2)2]V' 



(A5) 



which gives growth only when r = 2M + 0(e). For small 
initial values of (f> this means r ~ 2, but as </> grows 
because of the initial parametric resonance, the system 
goes out of resonance. 

We show that equations (A2,A3) can be solved exactly 
in terms of elliptic integrals. Multiply (A2) by a and 
(A3) by —6 and add to get: 



dt 



b^) = -{-)ab. 
r 



(A6) 



This equation is independent of the non-linear terms in 
(A2,A3); it would hold even if these terms were dropped. 
Note that exponential growth requires a, b to be of op- 
posite sign. The constraints expressed by equation (A6) 
allow the elimination of one degree of freedom: 



a = Acos'if, b = -Asin'if, 
with a relation between A and ^E": 

A = Aoexp / dt'(^)sin2*(t') 
Jo 2r 



(A8) 



with ^0 a-s an initial value. In the linear (Mathieu) case 
cos 5* is the constant — 4|/e, which yields the linear 
growth rate in equation (A5). But equations (A2,A3) 
yield two equations for the time evolution of A < ^P. 



The sum of these equations is a trivial identity, while the 
difference (using equations (A7,A8)) is: 



ecos2\I' - 2r* 



— -2-3Alexp 



)sin2*(t')- 
(A9) 



Now differentiate (A9) and use (A9) in the result to arrive 
at: 



2^- 



:{r'^ -4) sin2^'- 



2r2 



■ sin 45'. 



(AlO) 



This is readily checked to be an elliptic equation. We will 
not bother to study it here. All the physics is contained 
in the linearization of (AlO), which gives: 



*(t) = «'ocos[(A(t-to) 
with frequency 



A 



2r2 J • 



(All) 



(A12) 



Since — 4 is 0(e), so is A. The best case for growth is 

*(t) = (J)cosAt (A13) 

(so that A and b are equal initially). Evidently, from 
(A8) growth stops when ^E* = 0, or when t = 7r/2A. This 
means, as discussed in the main text, that growth cannot 
be unlimited. (However, when the frequency r grows sec- 
ularly, growth can continue unimpeded with a^,b^ ~ r, 
which maintains the resonant growth condition.) Gener- 
ally, no matter how small the initial values of the poten- 
tial (j), eventually becomes of order unity. The smaller 
the initial value, the longer this process takes. 
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